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ABSTRACT 
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tential. In particular we discuss possible signatures of halo wavefunctions 
in elastic and inelastic scattering experiments. 
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1 Introduction 



Reactions with radioactive beams are a useful tool to understand the properties of 
nuclei far from the stability line (see, e.g. ref. |]]). Up to the present these reactions 
are mainly restricted to studies of reaction cross sections and momentum distributions 
of the fragments @, ||. Elastic scattering has been studied in several experiments 
@) E| H 0' il • Inelastic scattering has also been studied, but very few experiments for 
Coulomb breakup processes || [K| only are available so far. Usual techniques with 
stable nuclei, like photo-nuclear absorption and electron scattering are far beyond the 
present experimental possibilities. The study of reaction cross sections, momentum 
distributions, and elastic scattering gives us a very limited access to the information 
on the internal structure of exotic nuclei, although many intriguing properties of these 
nuclei have been deduced in these experiments. An example is the discovery of the 
extended matter distribution in very light neutron-rich nuclei, so-called halo nuclei 

The study of inelastic excitation cross sections is the natural step to increase our 
knowledge on the nuclei far from the stability line. In fact, Coulomb and nuclear ex- 
citations in nucleus-nucleus scattering are well established tools for the spectroscopy 
of stable nuclei and are complementary to photo-nuclear and electron scattering ex- 
periments. Due to the low luminosity of radioactive nuclear beams, experiments are 
possible only when the cross sections are sufficiently large. For inelastic scattering 
this is the case for the Coulomb excitation of loosely-bound nuclei, e.g., 11 Li and 
11 Be, incident on heavy targets 10J. The beams of exotic nuclei are often available 
at intermediate and high energies, Ei a i, > 50 MeV/nucleon [JTJ . At these energies the 
Coulomb field favors the dipole excitations. On the other hand, the nuclear interac- 
tion favors the monopole and quadrupole excitations. 

A study of nuclear excitations of halo nucleus 11 Li has been done in ref. JTT]]. It 
was shown that the excitation of low-lying continuum states of u Li have substantial 
cross sections and it can reach some 100 mb/sr at forward angles. In this article 



we extend the model of ref. [11 to include Coulomb excitation processes. Inelastic 
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excitation processes are unavoidably accompanied by elastic scattering. We thus also 
make a study of elastic scattering processes with unstable nuclei. Under simplifying 
assumptions, which we try to justify in the next sections, elastic scattering allows one 
to get information on the ground state properties of the exotic nuclei, while the in- 
elastic scattering processes tells us about the transition densities for the excitation. If 
the ground state wave function can be deduced from the elastic scattering, the excited 
state wavefunctions are possible to infer from the inelastic excitation processes. 

The extraction of the above mentioned information from the experimental data 
suffers from several difficulties which arise due to the complexity of the reaction 
mechanism. The ideal situation occurs when the reaction mechanism can be explained 
with very simple models. In this article we use a simple and tractable model for the 
reaction mechanisms and apply it to the study of the scattering of unstable beams. 
Under some circumstances it is shown that our model can be very useful to extract the 
ground state properties of exotic nuclei. A study of this feature is presented in sections 
3 and 4. A brief explanation of the calculation ground state densities for exotic nuclei 
is done in section 3. These densities are used as inputs in the calculation of the elastic 
cross sections, which are compared to some available experimental data in section 4. 
To show the feasibility of the models used, we also compare the theoretical predictions 
with some experimental data for reactions with stable nuclei. 

The inelastic scattering cross sections are discussed in sections 5 and 6. In section 5 
we deduce formulas using the deformed potential model and folding model for inelastic 



excitations. The equations for the folding model has been used in ref. ||1 1|| . We also 
present appropriate formulas for the case of Coulomb excitation. These equations are 
applied in section 6 to discuss the angular distributions for inelastic processes with 
stable and unstable nuclear beams. In section 7 we present our conclusions. 
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2 Elastic Scattering 



The elastic scattering in nucleus-nucleus collisions is a well established tool for the 
investigation of ground state densities. This is because the optical potential can be 
related to the ground state densities by means of a folding of the nucleon-nucleon 
interaction with the nuclear densities of two colliding nuclei. But, this relationship 
is not quite straightforward. It depends on the effective interaction used, a proper 
treatment of polarization effects, and so on (for a review see, e.g., [0). At higher 
bombarding energies than Ex a & ~ 50 MeV/nucleon, a direct relationship between the 
nuclear densities and the optical potential is possible, as long as the effects of multiple 
nucleon-nucleon scattering can be neglected. The effects of real, or virtual, nuclear 
excitations are small since the excitation energies involved are much smaller than 



the bombarding energies. Lenzi, Vitturi and Zardi |L3| have performed an extensive 
study of the nuclear scattering of stable nuclei at high energies. From their study 
one concludes that a simple relationship between the ground state densities and the 
elastic scattering cross sections quite often yields very reasonable results, as compared 
to the experiments. We adopt here a similar approach and extend it to study the 
scattering of exotic nuclei. 

At high energies the elastic scattering cross section for proton-nucleus collisions is 



well described by means of the eikonal approximation [|14[]. The optical potential for 
proton-nucleus scattering is assumed to be of the form 



where 



and 



U{r) = U (r) + U s (r) (L.S) + U c (r) (1) 
U (r) = V R f R (r) - iW v / z (r) +4ia I W I ^fj(r) (2) 

Us(r) = 2 (^-) 2 V S - ±f 3 {r) (3) 
x -m 7T c / r ar 

are the central and spin-orbit part of the potential, respectively, and Uc(r) is the 
proton-nucleus Coulomb potential. The Fermi functions f\ are given by 

f i (r) = {l + ex 1 p[(r-R i )/a i ]}~ 1 , (4) 
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with Ri = riA 1/3 . 

In the eikonal approximation, the proton-nucleus elastic scattering cross section 
is given by |L4) 



da e i 



F{9) + G{6) 



where 



F(9) = f c (9)+ik J Q dbbJo(qb) exp[z X c(6) ^1-exp i Xo (b) 



cos 



(5) 



kb X s(b) (6) 



and 



G(9) = ik / dbb J\(qb) exp ixc(b) + iXo(b) 



sin 



kb Xs{b) 



(7) 



In the equation above q = 2&sin(0/2), where 9 is the scattering angle, Jo (Ji) is the 
zero (first) order Bessel function. The eikonal phase Xo(S) is given by 



1 r°° 

Xo(S)(b) = J U 0(s) (b, z) dz . 



(8) 



For the Coulomb eikonal phase we use the approximation, valid for a point nucleus 

2Z l Z 2 e 2 



Xc(b) = rj \n(kb) with r\ 



hv 



(9) 



where Z\ and Z 2 are the proton (Z\ = 1) and the nuclear charges, respectively. The 
Coulomb phase will be changed for a finite charge distribution of the nucleus ||15|| . For 
example, assuming a uniform charge distribution with radius R the Coulomb phase 
becomes 



Xc(b) 



rj \ Q(b- R) ln(fcfc) + Q(R - b) \n(kR) + ln(l + Jl- b 2 /R 2 ) (10) 



- Jl 



1 



byR 2 --(i-b 2 /R 2 f/ 2 
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where is the step function. This expression is finite for 6 = 0, contrary to eq. (9). 
If one assumes a gaussian distribution of charge with radius R, appropriate for light 
nuclei, the Coulomb phase becomes 



Xc (b)= V {Hkb)+ l -E 1 (b 2 /R 2 )} 



(11) 
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where the error function E\ is defined as 

Ei(x) = 



OO 

— dt . 

t 



(12) 



This phase also converges, as b — > 0. 

We have noticed in our numerical calculations that the above corrections to the 
approximation (9) do not modify the calculated cross sections appreciably. Moreover, 
one can show that in the eikonal approximation, the phase defined in eq. (9) yields 
exactly the Coulomb scattering amplitude 

Z x Z 2 e 2 



fc(0) 



exp | — it] In sin 2 (0/2) + + 2i0 o } 



(13) 



2[iv 2 sin 2 (9/2) 

where 0o — argT(l + ii]/2). This is convenient for the numerical calculations since eq. 
(6) is written with the separated contribution of the Coulomb scattering amplitude. 
Then, the remaining integral (the second term on the r.h.s. of eq. (6)) converges 
rapidly for the scattering at forward angles. If we use eqs. (10)and (11) the amplitude 
F(9) cannot be separated as in eq. (6) and the integral in h will converge very slowly. 

A more important correction, due to the Coulomb deflection of the trajectory, 
amounts to calculate all elastic and inelastic integrals (to be discussed in section 4) 
replacing the asymptotic impact parameter b by the distance of closest approach in 
Rutherford orbits, i.e., 

kb' = 7] + J v 2 + k 2 b 2 . (14) 



As shown by Vitturi and Zardi [IB] this correction leads to a considerable improvement 
of the eikonal amplitudes for the scattering of heavy systems. 

For nucleus-nucleus collisions the spin-orbit interaction, Us{r), and the surface- 
term of the imaginary potential (last term of eq. 2) are usually neglected. Whereas 
these terms are relevant for proton-nucleus scattering, they play no important role in 
nucleus-nucleus collisions. Thus, for nucleus-nucleus collisions, the scattering ampli- 
tude is given by 



dod 



F{9) 



where 



F{9) = f c (0) + ik | o °° db b J (qb) exp [i Xc (b)] {l - exp [%(&)] } 



(15) 



(16) 



Obviously, these equations can be obtained from eqs. (5-7) by setting xs — 0. 

A common way to relate the nuclear optical potential to the ground-state densities 
is to use the "tpp" approximation. This approximation has been extensively discussed 
in the literature | I2"| , |T3|j . In its simplest version, neglecting the spin-orbit and surface 
terms, the optical potential for proton-nucleus collisions is given by 

U (r) =< t pn > p n (r)+ < t pp > p p (r) (17) 

where p n (p p ) are the neutron (proton) ground state densities and < t pi > is the 
(isospin averaged) transition matrix element for nucleon-nucleon scattering at forward 
directions, 

1%%) 

t pi (q = 0) = -(2nh 2 /p) f pi (q = 0) = -— o pi (^ + i) (18) 

where a P i is the free proton-nucleon cross section and £ P j is the ratio between the 
imaginary and the real part of the proton-nucleon scattering amplitude. The basic 
assumption here is that the scattering is given solely in terms of the forward proton- 



nucleon scattering amplitude and the local one-body density [12 . 

For nucleus-nucleus collisions, we will use the same method which leads to an 
optical potential of the form 

U (R) = J < t NN (q = 0) > pi{R - r') p 2 (r') d 3 r' , (19) 

where R is the distance between the center-of-mass of the nuclei. We use the isospin 
average < t NN >= (t pp + t pn )/2. 

Without much computational effort the formula (19) is improved to account for the 
scattering angle dependence of the nucleon-nucleon amplitudes. A good parametriza- 
tion |HJ for the nucleon-nucleon scattering amplitude is given by /vAr(q) = (&aw /4vt) ctnn (i+ 
oinn) e~^ NNq2 . The nuclear scattering phase then becomes [jTJJ 

Xo (b) = J j dv dr' p x {v) lNN (\h - s - s'|) p 2 (r') (20) 

where the profile function 7atat(6) is defined in terms of the two-dimensional Fourier 
transform of the elementary scattering amplitude 

1nn(P) = — f exp [ - iq.bl /jv7v(q) rfq , (21) 
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and s, s' are the projections of the coordinate vectors r, r' of the nuclear densities 
on the plane perpendicular to the z-axis. For spherically symmetric ground-state 
densities eq. (20) simplifies to the expression 

Xn(P) = dqq p^q) f NN (q) p 2 (q) J (qb) , (22) 
Jo 

where Pi(q) are the Fourier transforms of the ground state densities. 

In the following we will use the eq. (17) in order to calculate the proton-nucleus 
eikonal phase-shifts. For nucleus- nucleus collisions we will use the formula (22), since 
we assume the spherical symmetry for the ground state distributions. We will now 
describe how we calculate the ground state densities for exotic nuclei. 



3 Ground State Densities 

We calculate the density distributions of halo nuclei based on the Hartree- Fock (H-F) 
approximation with Skyrme interaction. The Skyrme interaction is known to describe 
successfully the ground state properties (the binding energies, the rms radii and the 
charge distributions) of many nuclei in a broad region of the mass table |jTB| . 

It was pointed out that the separation energy of loosely bound neutrons plays an 



important role to study halo nuclei [|19], EH, ETJ. One needs an accuracy of about a few 



tens of keV to take into account the effect of the separation energy on the radii and 
the density distributions of halo nuclei, since they have extremely small separation 
energies as is shown in table 1. On the other hand, the H-F theory cannot provide the 
prediction for the separation energies within the accuracy of a few tens of keV |22] . It 
is known that higher order effects beyond the mean field approximation are necessary 



to describe more precisely the single-particle energies near the Fermi surface p3| . So 
far, several theoretical attempts [pJJ have been made in this direction, but the results 
do not satisfy the accuracy which is required for the study of halo nuclei. We take the 
following method [f20| , p5|j to improve the calculated separation energies rather than 
evaluating directly the higher order effects; the last neutron configuration is treated 
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differently from the other orbits in the H-F potential in order to reproduce properly 
the neutron separation energy of the nucleus. 

In the above procedure, the single-particle energy of the last orbit is adjusted to 
be the same as the empirical single-neutron separation energy in the case of odd-N 
system, while the empirical two-neutron separation energy is adopted for the last 
neutron orbit in the even-N nucleus. Although it is not obvious how much the two- 
body correlation between the halo neutrons affects the single-particle wave function, 
this choice for the even-N nucleus was pointed out reasonable to study the Coulomb 
dissociation cross sections n Li |Q. The soft dipole excitation in n Be JITJ is also 
described well by this wave function. 

The H-F equation for the Skyrme interaction can be written as 

r h 2 n 

-V— — V + V(r) Mr) = eMr) (23) 
l 2m* (r ) J 

where m*(r) is the effective mass. The potential V(r) has a central, a spin-orbit and 
a Coulomb term, 

V(T) ^central ^spin— orbit ^Coulomb- (24) 

This H-F potential can be expressed analytically in terms of the parameters of the 
Skyrme interaction [|18[. The central potential in eq. (2) is multiplied by a constant 
normalization factor f only for the last neutron configuration: 

Tr / \ „ Tr , N I f ^ 1 for last neutron con figuration /ri _ N 
V central (r) = fV H Mr), \ J f t j Qtherwise • (25) 



Numerical calculations are performed with the parameter set SGII[^] which gives sat- 
isfactory results for charge distributions of many nuclei, and also for the systematics 



of the nuclear radii in comparison with experimental data[27 



The empirical separation energies of light neutron-rich nuclei are tabulated in 
table 1. Nuclei with negative values in table 1 are unstable against neutron decay. 
We can see that the separation energies are extremely small in the cases of 6 He, 11 Li, 
11 Be and 14 Be which are known as halo nuclei. Calculated radii are also tabulated in 
table 2 together with empirical data. 
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4 Results for Elastic Scattering 



We apply the formalism presented in section 2 to the elastic scattering of stable 
nuclei in order to study the validity of the method. We will concentrate here in 
nuclei with spherical symmetry. The ground state densities of stable nuclei used here 
are parameterized as gaussian (G) and modified Fermi (MF) densities, as shown in 
table 3. The nucleon-nucleon cross sections used as input to construct the optical 
potentials in the "£pp" approximation are shown in table 3. A linear interpolation 
is done to find the appropriate set of parameters at energies between those shown 
in the table. For energies lower than 94 MeV/nucleon we use £nn = 0.5 fm 2 . As 
one can deduce from eq. (22) the g-dependence of f^N is very important with the 
parameters £,nn, as given in table 4 since the integrand in eq. (22) is only relevant 



for g < 1/6 l/vviVJV- The parameters for the densities of stable nuclei used are 
shown in table 3. 

In fig. 1 we show the experimental data from ref. for the elastic scattering of 
17 O + 208 Pb at 84 MeV/nucl. The curve is calculated using the nuclear phase-shift 
constructed as in eq. (22) and the scattering amplitude as in eq. (16). We observe 
that the agreement with the data is extremely good. The rainbow scattering at 
9 ~ 2.5° is also very well reproduced. However, one should be cautious with such an 
example since in this case the scattering amplitude is dominated by a sharp transition 
from no-absorption to strong absorption as the scattering angle increases. Basically, 
the optical potential has to have the feature of a large absorption at small distances 
and a small diffuseness. This behavior arises naturally within the framework of the 
"tpp" approximation for heavy systems. 

Light systems are more "transparent" and the collision at small distances (large 
scattering angle) are more sensitive to the details of the optical potential. To show 
this we plot in fig. 2 the elastic scattering of a + a at Ei a f, = 2.57 GeV as a function 
of the invariant momentum transfer t = —(pi +P2) 2 - The data are from ref. fl3l] . The 
calculated curve agrees reasonably well with the data at forward angles (small t) but 
deviates appreciably from it at larger angles. The forward scattering is dominated 
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by peripheral collisions for which multiple scattering is not relevant. Thus, we expect 
that the "tpp" approximation works well at forward angles. On the other hand, it 
has been shown that multiple collisions are very important for large scattering angles 



33|, so that the "tpp approximation fails to explain the data. 

We plot in fig. 3 the elastic scattering of 12 C+ 12 C at E iab = 85 MeV/nucl to study 
the model further. The data are from ref. [[34]]. The solid curve is obtained by using 
the eikonal approximation (eqs. 16 and 22) with a set of Woods-Saxon potentials 
which was constructed to give the smaller x-square fit to the data with a DWBA 
calculation using the code PTOLEMY EH]. This Woods-Saxon potential reproduces 



the data perfectly []34}| . The parameters obtained by this fit are 

V = -120 MeV, R v = 1.72 fm, and a v = 0.83 fm (26) 
W = -46.8 MeV, R w = 2.2 fm, and a v = 0.86 fm . 

The dashed curve is obtained with the "tpp" approximation which shows again the 
mismatch with the data at large scattering angles. We see that the "tpp approxima- 
tion gives a reasonable description of the elastic scattering only at forward angles. 
It is known that large scattering angles are affected by corrections due to real, and 



virtual (polarization) nuclear excitation p2 |. This means that a simple and unam- 



biguous relationship between the elastic scattering data and the nuclear ground-state 
densities does not exist at large scattering angles. We conclude that the disagreement 
between the experimental data and the method used here, for light systems and large 
scattering angles, is not a deficiency of the eikonal approximation, but of the con- 
struction of the optical potential (i.e. on the assumption of the validity of the "tpp" 
approximation) . 

The corrections due to multiple collisions using the Glauber multiple scattering 
series is rather cumbersome |3B| and a direct connection between the ground state 
densities of the nuclei and the scattering data is lost among a large variety of ap- 
proximations. We feel that, due to its simplicity, the "tpp" approximation is very 
appealing within its limitations. We therefore will use this approximation since the 
complications arising from a treatment of multiple scattering, and of other effects, 
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undermines our effort to extract information about the ground state densities of the 
nuclei from the elastic scattering. We will focus the following study mainly in the 
forward scattering region because of the limited validity of the model. 

Fig. 4(a) shows the elastic scattering of proton on 9 Li at Ei a t = 60 MeV. The 
data are from ref. [|J. The solid curve is obtained by using eqs. (5-7) with the 
parameters shown in table II of ref. ||] which was obtained by a %-square fitting 
to the data with a DWBA calculation. The dashed curve is obtained by using the 
optical potential constructed as in eqs. (17-19) together with the Hartree-Fock 9 Li 
ground-state density. In this case the surface and spin-orbit interaction are absent. 
The dashed curve clearly misses the experimental data. The absorption is greater 
than expected by the "tpp" method so that this difference cannot be ascribed only 
to the absence of the spin-orbit interaction. Quite a different scenario is presented in 
figure 4(b) where we plot the p + 11 Li scattering data at 62 MeV/nucleon from ref. 
M. The solid curve is calculated with the optical potential parameters of table II (set 
B) of ref. Q for eqs. (1-3), which were chosen so as to fit the experimental data. The 
dashed curve is obtained by using the eqs. (17-19) and the u Li Hartree-Fock ground 
state density, calculated as explained in section 2. The agreement with the data is 
quite good at forward angles, in contrast to the previous case. 

One might think that the good quantitative description obtained in figure (4b) is 
accidental, in view of the previous discussions. In order to clarify further the validity 
of the model, in fig. 5 we show the elastic scattering data of p+ 8 He at 72 MeV/nucleon 
together with a calculation using the optical potential constructed as in eqs. (17-19) 
by the solid curve. The ground state density of 8 He was calculated as explained in 
section 3. Although the data uncertainties at 9 > 40° do not allow for a test of 
the theory, the agreement with the data is almost perfect at forward angles. Such 
an agreement is very encouraging since the simple "tpp" approach is of very useful 
predictive power and it can be used to plan future experiments on radioactive beams 
scattering of protons. Also shown in fig. 5, the dashed curve is the scattering cross 
section obtained by using the 6 He ground state density. Two important features 
are seen. Firstly, the strength of the strong absorption in the p + 8 He system is 
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appreciably larger than in the p + 6 He system. Secondly, the diffraction minimum is 
at a larger angle for 6 He targets than for the 8 He ones, revealing a larger absorption 
radius of 8 He. 

The inclusion of the surface and spin-orbit terms of eq. (1) in a microscopic 
Glauber approach, as in eq. (17), to elastic scattering is rather complicated (see 
discussion in ref . JI^] ) . The absence of these terms in our approach may be considered 
as one of the main reasons for the discrepancies between the experimental data and 
the calculations. This problem should not occur for nucleus-nucleus collisions, since 
the folding of densities "smear out" the surface and spin-orbit terms of the nucleon- 
nucleus optical potential. 

Data on elastic scattering of halo nuclei has been taken during the last few years. 
Due to the poor energy resolution, the data are contaminated with inelastic scattering. 
In fig. 6, 7 and 8 we plot the data on quasi-elastic scattering of n Li+ 12 C at Ei ab = 637 
MeV, of 12 Be + 12 C at E lab = 679 MeV, and of u Be + 12 C at E lab = 796 MeV, 
respectively. The dashed curves are calculated by using the prescription given by eqs. 
(15-16, 22) and the Hartree-Fock densities calculated as explained in section 2. The 
data are from refs. [§, [/J. While the agreement is quite reasonable for the 12 Be + 12 C 
and 14 Be + 12 C data, it does not work so well for the 11 Li + 12 C data. The solid 
curves are calculated by including the inelastic excitation of the 2 + , and the 3~ states 
in 12 C which cannot be separated from the experimental data. In the next section 
we describe how to calculate the inelastic scattering excitation cross sections. We use 
the deformed potential model for the excitation of 12 C, with deformation parameters 
P 2 = 0.59 and j3 3 = 0.40 for the 2 + and the 3~ states, respectively. These are the 
same values used in refs. || [7|. We observe that, with the inclusion of inelastic 
scattering, the theoretical results for the 11 Li + 12 C collision is remarkably good, but 
the ones for l2 Be + 12 C and for 1A Be + 12 C are not, yielding too large cross sections 
at large angles. One possible reason for such a discrepancy is the the absence of 
dynamical polarization in our approach. The virtual excitations during the collision 
time is shown to affect the elastic scattering of halo neutron appreciably and has been 



studied by several authors [36]. An inclusion of such effects is beyond our model of 
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using a simplified approach to study the ground-state properties of the halo nuclei. 
Another possibility is that the folding prescription to determine the optical potential 
yields larger cross sections than the experimental ones, as we observed in fig. 3. 

We conclude that the simple "tpp" folding procedure is an useful method to de- 
scribe quantitatively the elastic scattering data of radioactive nuclei in most cases. We 
should also notice that the inelastic scattering makes it difficult to obtain good agree- 
ment with the data, especially at larger angles. However, a test of the ground-state 
densities of such nuclei is possible by comparing the data with the theoretical calcu- 
lations at forward angles. At large angles, especially with the presence of inelastic 
scattering, the results are however sometimes misleading. More detailed microscopic 
calculations would be necessary to obtain better agreement. In the opposite extreme, 
one may resort to the use of optical potential parameters to fit the data. In both 
cases, a link of the results to microscopic features of the nuclei is very difficult to 
achieve. 



5 Inelastic Scattering 
5.1 Nuclear Excitation 

We can use the Distorted Wave Approximation (DWBA) for the inelastic amplitude, 
assuming that a residual interaction U between the projectile and the target exists 
and is weak . The cross section for the excitation of a vibrational mode (A/i) with 
energy huj\ is given by 



where k\ is defined as 

*2,2 *2,2 

Ek x — — — — ^ = __2. _ %oj\ . (28) 

In eq. (28), M is equal to the reduced mass of the system and k is equal their relative 
momentum. The wavefunctions and \I/ describe the relative motion and the internal 
states, respectively. 
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In the particle-vibration-coupling model the transition matrix elements are given 

by 

M A „(r) =< * v |[/|* >= -y== Ux(r) (29) 

where 5\ = (3\R is the vibrational amplitude, or deformation length, R is the nuclear 
radius, and U\(r) is the transition potential. We will consider only low multipolarities, 
/ < 2 in the following. 

The deformation length 5\ can be directly related to the reduced matrix elements 
for electromagnetic transitions. Using well-known sum-rules for these matrix elements 
one finds a relation between the deformation length, and the nuclear sizes and the 
excitation energies. For isoscalar excitations one obtains 



Ti 2 < r 2 > ?ir h 2 1 

^k^r-^i ^ A < 2A+1 >ik (30) 



where A is the atomic number, < r 2 > is the r.m.s. radius of the nucleus, and E x is 
the excitation energy. 

The transition potentials for isoscalar excitations are 



U (r) = 3U opt (r) + r^^l , (31) 

for monopole, and 



= dUoM. , (32) 
dr 

for quadrupole modes. 

For dipole isovector excitations, the deformation length is given by 

r h 2 A 1 

TTipj NZ E x 

where Z (N) the charge (neutron) number. The transition potential in this case is 

„ l(r)=7 (^)(£k4*^), (34) 

where the factor 7 depends on the difference between the proton and the neutron 
matter radii as 

_ 2(iV — Z) _ R n — R p _ AR np 
3A 5 {Rn + Rp) Ro 
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Thus, the strength of isovector excitations increases with the difference between the 
neutron and the proton matter radii. This difference is accentuated for neutron-rich 
nuclei so that the isovector dipole excitations should be a good test for the quantity 
AR np which becomes very large for neutron-rich unstable nuclei. One can generalize 
this equation to higher isovector multipole excitations A > 2 by the substitution 

^ ^ Q<..» +Q « = Z(-i) A + [(l-i) A + (-l)^], (36, 

where Q^'^ are the effective charges (in units of e) of the neutron and the proton, 
respectively. 

An useful approximation, valid for hu\ <C E kx , is 



k ^ k «^-^w)= k » + -- (37) 
Further, by using the definition (27) and the eikonal approximation 

^^ ) ~exp{iq.r + j X (&)} , (38) 



we get 



M 



MO) = — ^ -== Sx i" Jn(2\ + 1) 



(A-/.)! 



(39) 



x 



poo poo / Z \ 

J o db b J^q t b) J_Jz P Am ( ^__ ) U x (b,z) , 

where and are the Bessel function and the Legendre polynomials, respectively, 
and q t = 2\ / k k\ sin(^/2). 

As seen in eq. (39), one needs to calculate two simple integrals to compute the 
inelastic scattering at intermediate energies with the deformed potential model. The 
scattering amplitudes will depend on the optical potential parameters and on the 
deformation length 5\. The deformed potential model is based on the assumption of 
a transition density peaked at the surface. Although this assumption is reasonable for 
the excitation of heavy nuclei ( e.g., 40 Ca, 208 Pb), it is rather crude for light nuclei, 
especially when the transition density extends radially beyond the nuclear size. This 
is the case for the soft multipole excitations, for which the transition densities have 
very long tails. 
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A more convenient way to describe the inelastic scattering of neutron-rich nuclei 
is to use the folding approximation. The assumption of a transition density peaked 
at the surface of the nucleus is not necessary. In this model the matrix element on 
the right-hand-side of eq. (27) is 

T x , = fd 3 Rf d 3 r <j>g*(R) U int (\R- r|) Sp x ,(r) 0«(R) , (40) 

where 5p XlI = ^* Xfl ^o is the transition density. U int (\R — r\j is the potential between 
each nucleon of the target and the projectile nucleus. Thus, the transition from 
the ground state to the excited state is directly calculated from the target nucleon- 
projectile interaction. 

For light targets, a gaussian parametrization of the (target nucleon)-projectile 
potential is adequate and yields simple formulas. This can be shown by means of the 
expansion 

U tnt (\R-r\) = (v + iw ) e -( R - r ) 2 A* 2 (41) 
= Mv + tw )e^ R2+ ^ 2 j A (2^) F V (R) r A ;(f), 

A/i ° 

where j\{ix) are the spherical Bessel functions calculated for imaginary arguments. 

Using this result in the eq. (27) and the definition Sp(r) = 8p x {r) Y x ^{i) we get 
(using R = (b, Z)) 



T Am = 4vr 3 / 2 (vo + iwo) 



A/i 

A/i 

POO 

X 



(2A + 1) (A -/x)! 



o 



\l (A + aO! 
db b J,{q t b) O x ,(b) 



(42) 



where 



POO 

O x »(b) = / dr r 2 5p x ,(r) F x ,(r,b) , (43) 

J 



with 



W ) =£^exp(-^) *,(f) (44, 
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A link between the deformed potential model and the folding model is obtained 
by using the approximation 



j 5 X dpo/dr /V2A+T, for A > 1; 
PX j So (3p + rdp /dr), for A = 0; [ > 

As in the deformed potential model, the scattering amplitude is determined by the 
optical potential parameters and the deformation length 5\. 

Another common approximation for 5p\ is provided by the Tassie model |5S| which 
gives 

^-VSTT(ir^< f° r A > 1 . (46) 
For A = 0, one uses eq. (45). In general, both models yield analogous transition 
densities for heavy nuclei and low collective states. 

These approximations assume that the transition density is peaked at the nuclear 
surface. As we will show later, this is a bad approximation for neutron rich nuclei. It 
is more convenient to use directly the transition density calculated from microscopic 
models and inserted into eqs. (42) and (43) to obtain the scattering amplitude. 



5.2 Coulomb excitation 



The subject of Coulomb excitation for heavy ion collisions at intermediate energies, 
has been discussed in ref . JJIJ , including effects of retardation and strong absorption 
effects. The Coulomb excitation amplitude for a given multipolarity EX is given by 



fEl,fi(E x , 



_ z V8n Z T eM aA , E^ 
3 h 2 ^hv a 

for El excitations, and 

fE2,n(E x ,9) = 



B(E1, E x ) 



1/2 f A ±1 (E X , 9), for/i = ±l; 
1 A {E X , 6), for fj, = 



2 0[ Z T eM aA , E^yi 

5 6 h 2 \hv a 



B(E2, E x 



1/2 



A± 2 (E X , 0)/ 7 , for n = ±2; 

x < -(2-v 2 /c 2 ) A ±1 (E X , 9), for/i = ±l; 
A (E X 9), for fj, = 



(47) 



(4? 
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for E2 excitations. 

In the above equations 7 = (1 — t^/c 2 )^ 1 / 2 , and /i is the azimuthal component of 
the transferred angular momentum. The functions A M are given by 

poo / E b \ v 1 

A,(E X , 6)= dbb J^qb) kJ-^-) exp \i x (b)] (49) 
Jo v 7ftf a / L J 

where q = 2fcsin(0/2), fc is the cm. momentum oi a + A, and J M (1^) are Bessel 
(modified) functions of order \i. x(b) = XN^b) + Xcip) is obtained from eqs. (8) and 
(9). 

Assuming that an isolated state is excited, and that it exhausts fully the sum 
rules, one gets (B(E\) = B{EX,E X )) 

B(E1) = ±jL^ e >, (50) 

v ; 4tt 2m N AE X ' v ; 

and 

„/™s fi 2 3i? 2 P f Z 2 /A, for isoscalar excitations: /r .^s 

B(E2) = e z x { . . , ., ,. (51) 

mjy 4:irE x Jvz/A, lor isovector excitations; 

The electromagnetic transition operators have a smooth dependence on the spatial 

coordinates (O(El) ~ r and 0(E2) ~ r 2 ). Thus, the electromagnetic excitation has 

a very weak dependence on the spatial form of the transition density and it does 

not seem to be a direct probe of the halo properties. However, the reduced matrix 

elements for halo nuclei are expected to be enhanced due to unique properties of the 

halo wavefunctions. For example, if the binding of the valence neutrons is very weak, 

as in the case of u Li, one expects a huge enhancement of the _B(_El)-strength at 

low energies (soft El-modes). This can be thought as due to the threshold effect of 

extended halo neutrons with respect to the core, which leads to large B(El)-values. 



6 Results for Inelastic 

We first apply the formulation of the 
of stable nuclei. In fig. 9 we plot the 



Scattering 

previous section to the inelastic scattering 
inelastic scattering of 84 MeV/nucleon 17 O 
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projectiles on lead. The excitation of the isoscalar giant quadrupole resonance in 
lead (E x = 10.9 MeV) is considered. Data points are from ref. [|3D]]. The curves 
are calculated by using the deformed potential model. The optical potential was 
calculated by using the folding procedure described in section 2. The deformation 
parameter @ 2 — 0.63 was used. This corresponds to 100% of exhaustion of the sum 
rule. The dashed curve is the contribution of the Coulomb excitation with B(E2) = 
0.73 e 2 b 2 , also corresponding to 100% of the sum rule (51). The dotted curve is 
the nuclear contribution only and the solid curve includes both contributions and 
their interference. We see that the scattering data are very reasonably described by 
this approach. The excitation of isovector electric dipole giant resonance in the same 
system is dominated by the Coulomb interaction and is also well described by eqs. 
(47,48), as shown in ref. []39|| . 

In general, both the deformed potential model and the folding model give very 
good agreement with the experimental data for the inelastic scattering of stable nuclei. 



This has been studied extensively, e.g., in ref. [37|. Also, in ref. [1 1] it was shown that 
a good agreement with the experimental data for the excitation of giant monopole 
and giant quadrupole resonances in lead by 172 MeV a's is obtained with the use of 
the eqs. (42-44), based on the eikonal approximation. The details of the oscillatory 
pattern of the angular distributions are easily understood in terms of the Bessel 
function J^(qtb) in eq. (42). 

An interesting application of this formalism is to the excitation of a Roper res- 
onance (E x = 500 MeV) in the a + p reaction with E a = 4.2 GeV. This is shown 
in fig. 10, together with the data points taken from ref. The data show a very 
steep angular dependence, characteristic of a monopole transition. As shown by those 
authors, the angular distribution can be well described by assuming that the Roper 
resonance is a monopole excitation, exhausting a large fraction of the monopole sum 
rule, eq. (30). The solid curve in the figure is calculated by using a folding po- 
tential derived according to eq (17) and the transition amplitude calculated within 
the deformed potential model, eq. (39), and using 5 calculated as in eq. (30), with 
< r 2 >= 3R 2 /5 = 0.69 fm 2 (and A = 1). The agreement with the data is remarkably 
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good and the value of < r 2 > is in reasonable agreement with the predictions for the 
electromagnetic radius of the nucleon in models which include sea-quark polarization 



effects 

To our knowledge, up to the present date very few experimental data exist on the 
inelastic scattering of exotic nuclei. Using the formalism in eqs. (40-44) predictions 
for the inelastic excitation of soft monopole and quadrupole modes in the 11 Li incident 
on carbon targets has been done |T|. It was shown that the particular oscillatory 
pattern of the excitation cross sections could be an important tool to identify different 
excited states in these nuclei. This is a well known method in reactions with stable 
nuclei, which allows one to distinguish, e.g., monopole from quadrupole excitations. 
The extension of the halo in loosely-bound nuclei would also affect appreciably the 
drop in magnitude of these cross sections with increasing scattering angle . 



In figure 11 we plot the excitation cross section for the 2 + state (E x = 3.57 MeV) 
in 8 He by 72 MeV protons. The data are from ref. ||. We used the folding model 
for the calculation of the optical potential. The inelastic cross section was calculated 
with the deformed potential model and a deformation parameter = 0.32. As seen, 
the agreement with the data is quite good. 

Finally we discuss the Coulomb excitation of exotic nuclei. As an example we 
consider the Coulomb excitation of 11 Be projectiles with 45 MeV/nucleon incident 
on lead targets. We consider the electric dipole excitation of the l/2 + ground state to 
the l/2~ state at 0.32 MeV. This experiment has been recently done at GANIL j|2 



as an initiative of studying the properties of low-lying states in exotic nuclei. We use 
eq. (47) with B(E1) = 0.116 e 2 fm 2 which is deduced from the lifetime measurement 
by Millener et al. P3| . We obtain the solid curve shown in fig. 12. The Coulomb 
cross section is peaked at very forward angles, as expected. The angular integrated 
cross section is equal to 210 mb. We omit the nuclear contribution which is very small 
and only important for large scattering angles. The cross section presents a wiggling 
at large scattering angles, characteristic of diffraction patterns in inelastic scattering. 
It is instructive to compare this result with a semiclassical calculations, based on 
Rutherford orbits for the trajectory of the nuclei and time-dependent perturbation 
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theory In this case, the excitation cross section is given by 



da 
dtt 



16 

9hc 



7T 3 B(E1,E X ) 



driEi 



(52) 



where the equivalent photon numbers driEi/dQ are given analytically by 



dn E i 



dn, 




T (-) 2 ^ C 2 e"< { 




[K K {<)f + K C «)] 2 } (53) 



where e = 1/ sin(6 l /2), a = 1/137, £ = E x ao/ , yhv, with a = Z a ZAe 2 /2EL a b- 

Using the above expression for dnEi/dQ a we obtain the dashed curve in figure 
12 for the same reaction One observes that the quantum result does not deviate 
from the semiclassical result appreciably The good agreement especially at lower 
angles is quite satisfactory. This shows that strong absorption is not relevant for the 
scattering at low angles. The peak at small angles is a consequence of the adiabaticity 
condition. For ( ^> 1 {6 a <C 0.2°) the Coulomb field is too weak to provide the 
necessary excitation energy. For ( < 1 (9„ » 0.2°) the Coulomb field is too strong 
and privileges the excitation of the projectile to higher energy states. Therefore, the 
cross section at a fixed relative energy of the fragments in the final channel has a peak 
at the optimal scattering angle corresponding to that energy. For the case above this 
angle is about 0.2°. 

7 Conclusions 

We have studied the applicability of simple concepts from scattering theory for inter- 
mediate energy collisions as a tool to study the ground state densities and the tran- 
sition probabilities in reactions with radioactive nuclei. The eikonal approximation 
together with the "tpp" approximation yields the simple and transparent formulas. 
This model gives very reasonable results for elastic scattering cross sections at forward 
angles which can be used to test the ground state densities of radioactive nuclei. The 
extended nuclear matter in exotic nuclei is manifest in the magnitude of the elastic 
cross sections as well as in the position of the first minimum. 
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The transition densities from the ground state to continuum states have the char- 
acteristic of an extended tail at large distances from the nuclear center. As discussed 



in ref. [JLlj] the folding model for the inelastic scattering cross sections are more appro- 
priate for the study of reactions with exotic nuclei. In this article we have presented 
a few more cases, complimentary to those studied in ref. [|ll|]. Comparisons with 
the experimental data in figs. (9-12) shows that the formalism presented here is well 
suited to the study of inelastic excitation of radioactive beams. The cross section 
for some cases predicted here are large and are well under the experimental possibil- 
ities. Of special interest would be a possible identification of soft multipole modes in 
such experiments. This could be accomplished by looking at the angular pattern of 
inelastic scattering since different multipolarities yield quite different patterns. 
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Figure Captions 



Fig. 1 Elastic scattering cross section of 17 O + 208 Pb at 84 MeV/nucl. The data 
points are from ref. P0| . The solid curve is a calculation with the "tpp" approxima- 
tion. 

Fig. 2 Elastic scattering cross section of a + a at Ei a b = 2.57 GeV as a function 
of the invariant momentum transfer t = — (pi + P2) 2 - The data points are from ref. 



| 3jJ . The solid curve is a calculation with the "tpp" approximation. 



Fig. 3 Elastic scattering cross section of 12 C + 12 C at .E^ = 85 MeV/nucl. The 
solid curve uses the eikonal approximation and the WS optical potential parameters 
given by eq. (26). The dashed curve uses the "tpp" approximation. The data are 
taken from ref. [04i . 



Fig. 4 (a) p + 9 Li elastic scattering at E p = 60 MeV. Data are from ref. Q . (b) 
p + 11 Li elastic scattering at E p = 62 MeV. Data are from ref. Bj]. 

The solid curves are obtained by using the optical potential as in eqs. (1-3) with 
the same set of parameters used in ref. £§]. The dashed curves use "tpp" potentials 
constructed as in eqs. (17-19). 

Fig. 5 Elastic scattering cross section for p + 8 He at 72 MeV/nucl. The solid 
curve is a calculation using eqs. (17-19), while the dashed curve is for the system 
p + 6 He at the same bombarding energy. The data are from ref. || . 

Fig. 6 Quasi-elastic scattering of ll Li + 12 C at Ei a b = 637 MeV. Data are from 
ref. ||. The inelastic contribution to the cross section was added to the elastic cross 
section (solid figure). Pure elastic cross section is given by the dashed curve. The 
deformation parameters fti = 0.59 and fa = 0.40 for the 2 + and 3~ states in 12 C were 
used. 

Fig. 7 Quasi-elastic scattering of l2 Be+ 12 C at Ei^ = 796 MeV. Data are from 
ref. 0. The inelastic contribution to the cross section was added to the elastic cross 
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section (solid figure). Pure elastic cross section is given by the dashed curve. The 
deformation parameters fa = 0.59 and (3^ = 0.40 for the 2 + and 3~ states in 12 C were 
used. 

Fig. 8 Quasi-elastic scattering of 14Be+ 12 C at Ei a t = 679 MeV. Data are from 
ref. 0. The inelastic contribution to the cross section was added to the elastic cross 
section (solid figure). Pure elastic cross section is given by the dashed curve. The 
deformation parameters fa = 0.59 and (3^ = 0.40 for the 2 + and 3~ states in l2 C were 
used. 

Fig. 9 Excitation cross section of the giant quadrupole resonance in 208 Pb in 
the reaction 17 + 208 Pb at 84 MeV/nucleon. Data are from ref. [pOj] . The optical 
potential was calculated by the "tpp" approximation. The deformed potential model 
is used to calculate the nuclear excitation cross section which is drawn by a dotted 
curve. The dashed curve is the Coulomb excitation cross section. The total cross 
section due to Coulomb and to nuclear excitations is shown by the solid curve. 

Fig. 10 Angular dependence of the inelastic cross section for the excitation of 
the Roper resonance in proton with the reaction a + p at E a = 4.2 GeV. Data points 



are from ref. |40l . The solid curve is calculated by using the deformed potential 
model. 

Fig. 11 Excitation cross section of the 2 + -state (E x = 3.57 MeV) in 8 He by 
protons with 72 MeV. Data points are from ref. 0. The solid curve is obtained 
by using the "tpp" potential for p + 8 He and the deformed potential model for the 
excitation, with 02 = 0.32. 

Fig. 12 Coulomb excitation cross section of the l/2~ state at 320 keV in 11 Be 
incident at 45 MeV/nucl. on a lead target. The solid curve uses the eqs. (47-49), 
with the eikonal phase calculated with the "tpp approximation. The dashed curve is 
a semiclassical calculation using eqs. (52-53). 
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Table 1. Single-neutron and two-neutron separation energies of neutron-rich nuclei. 
Data are taken from ref. |28| . 



A 

A 


T7T 

J 


S n (MeVj 


O /TV K A A 

5 2 n(MeV) 


6 He 


0+ 


1.86 


0.97 


8 He 


0+ 


2.58± 0.01 


2.14± 0.05 


9 Li 


3~ 
2 


4.06 


6.10 


n Li 


3~ 
2 


0.73± 0.05 


0.31± 0.05 


n Be 


1 + 
2 


0.51 


7.32 


12 Be 


0+ 


3.17 


3.67 


14 Be 


0+ 


3.35± 0.11 


1.34± 0.11 



Table 2. Calculated mass radii and observed interaction radii of halo nuclei. Radii of 
some halo orbits are also tabulated in the table. The core radii are obtained by the 
H-F calculations, while the values for the halo configurations are calculated by using 
the renormalized potential (3). Data are taken from refs. ||. 



A 


ilast 


s/(r 2 )aai (fm) 


\/i r2 )ex P (fm) 


6 He 




2.66 


2.48± 0.03 


8 He 




2.59 


2.49±0.03 


9 Li 




2.45 


2.41± 0.02 


n Li 


lPl/2 


5.36 








3.08 


3.20± 0.03 


n Be 


2si /2 


6.29 








3.01 


2.86± 0.04 


12 Be 


2si /2 


3.46 






2.68 


2.82± 0.04 


14 Be 


2si/ 2 


5.28 






3.16 


3.33± 0.17 
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Table 3. Parameters of the ground states densities of stable nuclei which are given 
by either gaussian (G) or modified fermi (MF) distributions. Data are taken from 



Nucleus 


Model 


R 


a 


c 


a 


Gaussian 


1.37 






n c 


MF 


2.34 


0.5224 


-0.149 


120 


MF 


2.34 


0.5224 


-0.149 




MF 


2.61 


0.513 


-0.051 


Vl O 


MF 


2.61 


0.513 


-0.051 


208 


MF 


6.62 


0.549 






Table 4. Parameters [17] for the nucleon-nucleon amplitude, fNN(q) = {J^nn / 4tt ) ct^n 
a NN ) e 



-S.NNQ 



E [MeV/nucl] 


a NN [fm 2 ] 


&NN 


£ NN [fm 2 ] 


30 


19.6 


0.87 




38 


14.6 


0.89 




40 


13.5 


0.9 




49 


10.4 


0.94 




85 


6.1 


1 




94 


5.5 


1.07 


0.51 


120 


4.5 


0.7 


0.58 


200 


3.2 


0.6 


0.62 


342.5 


2.84 


0.26 


0.31 


425 


3.2 


0.36 


0.24 


550 


3.62 


0.04 


0.062 


650 


4.0 


-0.095 


0.08 


800 


4.26 


-0.075 


0.105 


1000 


4.32 


-0.275 


0.105 


2200 


4.33 


-0.33 


0.13 



29 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502019vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502019vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502019vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502019vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502019vl 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502019vl 



